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Abstract: We apply the time-renormalization group approach to study the effect of pri- 
mordial non-Gaussianities in the non-linear evolution of cosmological dark matter density 
perturbations. This method improves the standard perturbation approach by solving renor- 
malization group-like equations governing the dynamics of gravitational instability. The 
primordial bispectra constructed from the dark matter density contrast and the velocity 
fields represent initial conditions for the renormalization group flow. We consider local, 
equilateral and folded shapes for the initial non-Gaussianity and analyze as well the case 
in which the non-linear parameter /nl parametrizing the strength of the non-Gaussianity 
depends on the momenta in Fourier space through a power-law relation, the so-called run- 
ning non-Gaussianity. For the local model of non-Gaussianity we compare our findings 
for the power-spectrum with those of recent N-body simulations and find that they accu- 
rately fit the N-body data up to wave-numbers k ~ 0.25/i/Mpc at z = 0. We also present 
predictions for the (reduced) matter bispectra for the various shapes of non-Gaussianity. 
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1. Introduction 

During the last years, the use of semi-analytic methods to study the Large Scale Structure 
(LSS) formation via gravitational instability has experienced a renewed interest, moti- 
vated mainly by high precision measurements of statistical properties expected in the next 
generation of galaxy surveys. These methods typically involve resummation or closure 
prescriptions and need to be compared with numerical simulations to assess their range 
of validity. Renomalized perturbation theory |Q, |2| is based on suitable resummations of 
perturbative terms and is expected to work well for the linear and "weakly non-linear" 
regime. Perturbation theory can be rewritten in a compact way which makes it possible 
the use of standard tools of quantum field theory, thus allowing a systematic treatment of 
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perturbative terms in terms only of a finite set of basic building blocks, namely, the initial 
fields, the linear propagator, which describes the linear evolution of the initial fields and 
an interaction vertex, responsible for the non linear mode coupling of the fields. This leads 
to the resummation of infinite classes of perturbation theory diagrams. The renormaliza- 
tion group perturbation theory |^] attempts to regulate the relative divergence of one-loop 
standard perturbation theory (for a review, see using renormalization group methods. 
A path-integral formulation of the Vlasov equation has been developed in in terms of 
the distribution function in the phase space. The lagrangian resummation theory Q is an 
extension of the well-developed Lagrangian perturbation theory. Refs. 0] and ||8| proposed 
a formalism based on path-integrals and renormalization group techniques which provides 
a systematic resummation scheme of cosmological perturbation theory (subsequently, this 
formalism was extended in Ref. ^ to calculate the non-linear propagator in the presence of 
non-Gaussian initial conditions of the fields). Finally, in the Time- Renormalization Group 
(TRG) approach proposed in Ref. the power spectrum, the bispectrum and higher 
order correlations, are obtained - at any redshift and for any momentum scale - by inte- 
grating a system of differential equations. The method is similar to the familiar BBGKY 
hierarchy. Truncating at the level of the trispectrum, the solution of the equations corre- 
sponds to the summation of an infinite class of perturbative corrections. The approach can 
be seen as a particular formulation of the renormalization group, in which time is the flow 
parameter. Compared to other resummation frameworks, this scheme is particularly suited 
to cosmologies other than ACDM and has been recently applied to compute the non-linear 
spectrum in presence of neutrino masses A critical look at the various cosmological 
perturbation theory techniques can be found in . 

On the other hand, over the last decade a great deal of evidence has been accumulated 
from the Cosmic Microwave Background (CMB) anisotropy and Large Scale Structure 
(LSS) spectra that the observed structures originated from seed fluctuations generated 
during a primordial stage of inflation. While standard single-field models of slow-roll infla- 
tion predict that these fluctuations are very close to Gaussian (see Q), non-standard 
scenarios allow for a larger level of non-Gaussianity (NG) (see and references therein). 
A signal is gaussian if the information it carries is completely encoded in the two-point cor- 
relation function, all higher connected correlators being zero. Deviations from Gaussianity 
are therefore encoded, e.g., in the connected three- and four-point correlation functions 
which are dubbed the bispectrum and the trispectrum, respectively. A phenomenological 
way of parametrizing the level of NG is to expand the fully non-linear primordial Bardeen 
gravitational potential <I> in powers of the linear gravitational potential <I>l 



The dimensionless quantity /nl sets the magnitude of the three-point correlation function 
1 15]. If the process generating the primordial NG is local in space, the parameter /nl 
in Fourier space is independent of the momenta entering the corresponding correlation 
functions; if instead the process which generates the primordial cosmological perturbations 
is non-local in configuration space, like in models of inflation with non-canonical kinetic 
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terms, /nl acquires a dependence on the momenta. It is clear that detecting a significant 
amount of NG and its shape either from the CMB or from the LSS offers the possibility 
of opening a window into the dynamics of the universe during the very first stages of its 
evolution. Non-Gaussianities are particularly relevant in the high-mass end of the power 
spectrum of perturbations, i.e. on the scale of galaxy clusters, since the effect of non- 
Gaussian fluctuations becomes especially visible on the tail of the probability distribution. 
As a result, both the abundance and the clustering properties of very massive halos are 
sensitive probes of primordial non- Gaussianities H, |l7|, IH, H, 10, 13, H, H, 11, H, and 



could be detected or significantly constrained by the various planned large-scale galaxy 
surveys, both ground based (such as DES, PanSTARRS and LSST) and on satellite (such as 



EUCLID and ADEPT) see, e.g. [||] and |2§. Furthermore, the primordial NG alters 



the clustering of dark matter (DM) halos inducing a scale-dependent bias on large scales 



|2S, 3C, 31, |3^ while even for small primordial NG the evolution of perturbations on super- 



Hubble scales yields extra contributions on smaller scales 1 33 , l2H| . The strongest current 



limits on the strength of local NG set the /nl parameter to be in the range —4 < /nl < 80 
at 95% confidence level 



In this paper, we implement the TRG approach to investigate the effects imprinted by a 
primordial NG in the non-linear evolution of cosmological perturbations. In particular, we 
will compute the DM power spectrum and the (reduced) bispectrum when some primordial 
initial NG condition is present. In the TRG approach this information is promptly encoded 
in the RG equations for the power spectrum and the bispectrum. As a consequence, we 
can easily study the impact of the various shapes of NG. In particular, we consider the 
local, the equilateral and the folded shapes. As we shall see, our approach is based on 
a closure assumption, i.e. we solve the RG equations for the power spectrum and the 
bispectrum, setting the connected n-point correlators to zero starting from the four-point 
correlator, the trispectrum. While this assumption represents a clear limitation of our 
approach, nevetheless it provides a well-defined and controllable scheme which allows us to 
quantify in a precise form the theoretical error associated with the method. For the local 
model of non-Gaussianity we will compare our findings for the power-spectrum with those 
of recent N-body simulations and find that they accurately fit rather well the N-body data 
up to wave-numbers k ~ 0.25/i/Mpc. We will also present predictions for the (reduced) 
matter bispectra for the various shapes of NG. 

The paper is organized as follows. In Section 2 we review the time evolution of the 
correlators governed by Eulerian dynamics and discuss our truncation hypothesis in com- 
parison with other approaches in the literature. In Section 3 we consider the primordial 
non-Gaussianities compatible with the truncation scheme of this paper and some theo- 
retically motivated models for the primordial bispectrum. In Section 4, we introduce the 
reduced bispectrum as a relevant tool to study the effects related to the shape of the pri- 
mordial bispectrum. In Section 5 we present our results for the power-spectrum and a 
comparison with the results of N-body simulations for "local" NG. We also present several 
illustrative plots of the reduced bispectrum for each model of NG. In Section 6, we end 
with our conclusions and discuss future perspectives on this line of work. 
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2. Dynamics of Gravitational Instability 



2.1 Eulerian Dynamics 

Our starting point are the hydrodynamic equations in the "single stream" approximation 
of a self-gravitating fluid made of cold dark matter collisionless particles in an expanding 
Universe. In terms of the mass-density fluctuation 5, the peculiar velocity v of the fluid 
and the peculiar gravitational potential (f), the dynamics of the system is governed by the 
system of equations 

|^ + V-[(l + <5)v] = 0, (2.1) 
^ + + (v V)v = - V0, (2.2) 

= ^n^a^6. (2.3) 

The first two equations are the continuity and Euler equations respectively, while the third 
one is the Poisson equation obeyed by perturbations on sub-horizon scales. Here, a is the 
scale factor of the background, t = J dt/a is the conformal time, H = dloga/dr = aH 
and is the matter density parameter. In the following we restrict our discussion to an 
Einstein-de Sitter model, so we take 0^, = 1^. As usual, we take the divergence of the 
Euler equation and define the velocity divergence ^ = V • v and Fourier transform. The 
resulting equations can be written in a compact notation as 

{^abdri + ^ab)Vb{K V) = e'^labciK "P, -q)'^b(P, V)'Pc{(l, V), (2-4) 

where we define the two-component field ipa{a = 1,2) through 

(::i:;:!)--<-'(-.t?/«)^ 

= a/ttin for an initial scale factor am conveniently fixed at an early epoch, 

"=(-3/2372) P-") 

and 7abc is a vertex function which describes the non-linear coupling of the modes and 
whose only non-vanishing components are 

7i2i(p, q, k) = h{p + q + k)^^^:^, (2.7) 

7222(P, q, k) = i<5(p + q + k)^5L±^^qJ^ , (2.8) 
7ii2(p, q, k) = 7121 (p, k, q). (2.9) 
Summation (integration) over repeated index (momenta) is understood in Eq. (|2.^). 



^For a generalization to more general cosmological backgrounds the reader is referred to Ref. |ic| ]. 
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2.2 Time evolution of the correlators 

Following Ref. ]1C1U the time evolution of the field correlators can be obtained directly from 
iterative application of the equations of motion ( |2.4[) . The result is an infinite tower of 
coupled differential equations relating the evolution of n-point correlators with the n and 
(n + l)-point correlators evaluated at the same time: 

+e'^Jacd{V^cfdfb) + e'^Jbcdifafcfd) , 

dr, ifaVb^c) = -^adiVdVb^c) " ^bdi^afdfc) " ^cdifafbfd) 
+e'^lade{VdVeVbVc) + e'^-fbdeiv^a^d^e^c) 
+e'^Jcde{^a'PbVdVe) , 

dr, ifafb^cfd) = ■■■ 

(2.10) 

Here, we omitted the time and momentum dependence in order to have compact expres- 
sions. To solve this system of equations is equivalent to apply the standard perturbation 
theory approach in which one calculates the n-point correlation functions by summing an 
infinite series of perturbative corrections depending on the interaction vertex 7, the red- 
shift and the initial statistic of the fields. Next, we introduce the following nomenclature 
for the first correlation functions of the fields 

((y9a(k,7?)(/?b(q,r/)) = 5z)(k + q)Pafe(k,r/) 

{ipa{Kr])^b{<l,v)^c{p,v)) = ^Z)(k + q+ p)Sa6c(k, q, p; v) 

(V3a(k, r?)<^6(q, v)'Pcip, v)'Pd{r, r/)) =6d{^ + q)(^D(p + r)P„6(k, r?)Pcd(p, 

+5i5(k + p)<5D(q + r)Pac{K v)Pbd{cL, v) + Sd{^ + r)5D(q + p)PadiK v)Pbc{<i, v) 

+Sd{^ + q P + r)Tabcd{K q, P, r; r/), 

where Pab is the power spectrum (PS) , Babe is the bispectrum (BS), and Tabcd is the 
connected part of the four-point function, the trispectrum. As said before, if we want to 
calculate the time evolution of the power-spectrum it will involve the time evolution of the 
bispectrum, which in turn involves the trispectrum and so on, giving us the infinite tower 
of equations ( p.lOj) . This system can be truncated if we neglect the trispectrum Tabcd = 0, 
letting us with the following system of two equations for the PS and the BS: 

dr, PabO^ , V) = -^^ac(k , ri)Pcbi'k , Tj) - Qbci^ , ??)-Pac(k , T]) 

+e'' j d^q [7acd(k, -q, q - k) Bbcdi}^, -q, q - k; ??) 

+Bacd{K -q, q - k; r?) 7f,cd(k, -q, q - k)] , 
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dr, BabciK -q, q - k; r?) = -Oad(k , r])Bdbc{K -q, q - k; ry) 

-^bdi-fi,r])Badc{K -q, q-k; ??) 
-Ocd(q - k , r])Babd{K -q, q - k; r/) 
+2e'' [7ade(k, -q, q - k)Prf6(q , rj)Pec(k - q , ??) 
+76de(-q, q-k, k)Pdc(k - q,r/)Pea(k,r?) 
+ 7cde(q-k, k, -q)Prfa(k,77)Peb(q,??)] • (2.11) 

The system can be formally solved in the form 

-Pafe(k , r/) = Qaciv, ^)gbd{ri, 0)Pcd(k , ?/ = 0) 

+ y drj'e'^' j d^qgae{v,v')9bf{v,v') 

X [7ecd(k, -q, q-k) B/cd(k, -q, q-k; ry') + (e ^ — ^ /)] 
BabciK -q, q-k; r/) = gad{v,0)9be{r],0)gcf{r],0)Bdef{K -q, q-k;r/ = 0) 

+2 / drj'e"^ gad{v,v')9be{v,v')9cf{v,v') 
Jo 

X [7dgh(k, -q, q - k)Peg(q,??')-P/fe(k - q,?/') 

+7egfc(-q, q-k, k)Pfg{k - q,r]')Pdh{'k ,r]') 

+ jfgh{ci-k, k, -q)Prfg(k,r?')Peh(q,??')] , (2-12) 

where gac{k;r],r]') is the linear propagator which is the Green's function of the linearized 
version of equation ( |2.4| ) (by setting 7 = 0) and gives the time evolution of the fields at 
linear order: 93^(k, r/) = gac{k.,ri,ri')(p^{'k,r]'). The subscript L stands for the linear order 
approximation. From Eq. ( p.l2D one gets an insight on the effect of our only approximation, 
namely, the truncation prescription Tabcd = 0. First, the solution for the bispectrum is 
formally "tree- level" , i.e. it contains no momentum integration. It has the same structure 
as the lowest perturbative contribution to the bispectrum, but with the linear power spectra 
replaced by the resummed, time-dependent ones. Inserting the bispectrum in the solution 
for the power spectrum, we see that it is formally one loop (one momentum integral), but 
with resummed power spectra in the loops instead of linear ones. In diagrammatic terms, 
starting from the one loop diagrams for the power spectrum (which have been computed by 



Taruya et al. [35|, for the non-Gaussian case), one is adding all the infinite contributions 
which can be obtained by adding corrections to the internal power spectra lines at any 
perturbative order. We use the language of Feynman diagrams introduced in [0, ^ |l^ to 
represent diagrammatically the perturbative contributions to the correlators. The basic 
building blocks of this diagrammatic approach are shown in figure [l|. In terms of them, 
we can represent the lowest order terms of this procedure as shown in figure |2| On the 
other hand, the comparison of the bispectrum with the computation by Sefusatti [^6| is not 
so straightforward. Compared to that computation we, again, add the corrections to the 
power spectra lines at all orders, but we do not include the corrections to the vertex, which 
Sefusatti computed at one-loop. In order to do that, we should go beyond our truncation 
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approximation and include the trispectrum in the hierarchy of our equations. This can 
be easily understood if one recalls that in perturbation theory the connected four-point 
correlator always gets a contribution which is generated by the bispectrum. It is precisely 
this contribution which generates the one-loop correction to the bispectrum which can be 
seen as a correction to the interaction vertices 7. Diagrammatically, the terms included 
and not included in this approximation are shown in figure ^. Including the trispectrum 
is straightforward in the TRG approach and we leave this to future work. On general 
grounds, we expect that the inclusion of the trispectrum will push the length scale where 
the method works to smaller scales. For a more detailed explanation of the diagrammatic 
of the method the reader is referred to Sections 3 and 4 in Ref. [|lO| ]. 




e'^ yabc Bkbc 



Figure 1: Basic building blocks for the perturbative approach with non-gaussian initial conditions. 
The linear propagator gabiVa,Vb) is causally oriented from the time r]b (dashed segment) to rja (solid 
segment). The interaction vertex e^jabc introduces the non linear mode coupling of the fields. 
The square on the linear power spectrum diagram, represents the initial conditions of the power 
spectrum, while the triangle in the linear bispectrum diagram represents the initial bispectrum, 
given by the sum of the newtonian contribution at z = Zi„ and the contribution due to primordial 
non-Gaussianities. 



3. Bispectrum Initial Conditions 



In this section we ca/lculate the initia.1 conditions, at redshift z — '^^m; 

for the system ( |2.11] ) 



for our case of non-Gaussian primordial statistics. To this end, we calculate the three-point 
correlators for the matter density contrast and velocity divergence up to second order terms 
taking into account the effect imprinted by primordial non-Gaussianities. In other words, 
we assume that second order perturbation theory holds from z = 00 to z = Zin, and the 
resummation embodied by the TRG is active from Zm down to 2; = 0. We follow the 



notations and results of Ref. |33| where the effects of primordial non-Gaussianities in the 



general relativistic cosmological evolution of matter perturbations were calculated up to 
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--0- - - -■■ 




Pab 



+ 3 



-<;:>- 



Babe 




+ cyclic permutations 



Figure 2: Diagrammatic representation of the formal solution for PS and BS. The first two lines 
represents the lowest order contributions to the PS equation. The first line is the 0(7^) contribution 
to the PS. The second line is the contribution coming from the non- vanishing initial BS which is a 
tree level, 0(7), contribution. The the third line is the tree level, 0(7), contribution to the BS. 




yes 




no 



Figure 3: Pcrturbative corrections to 0(7'^) order for the bispectrum which are included (left), 
and not included (right) by implementing the truncation prescription Tabcd = 0. 



second order in both the comoving and Poisson gauge, deep in the matter dominated epoch. 



Under certain conditions that we shall specify shortly, the results of [33] can be ex- 
tended to the case of a ACDM cosmology. In this case, the linear order solution of the 



fields are S^^^t) = D+{a)6^^' and el'>{T) = -Uf{T)D+{a)5)^' where D+{a) is the linear 
growth factor of density perturbations and /(r) = dlnL'+(a)/dlna. Information about 
the cosmological parameters of the background will be encoded in -C)+(a) and /(r). At this 
point, being interested in fixing the initial conditions at some early we use the approx- 
imation rim/ /{t)"^ ~ 1 which, as discussed in Ref. [@](see also Appendix B of Ref. [0), is 
a condition respected during most of the cosmological time evolution. Within this approx- 
imation one can simply use the kernels as can be derived from Ref. , properly replacing 
the conformal time r in terms of the functions D^{a) and /(r). Then, the second order 



.(0) 



3(1) 



(0) 
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expansions for the matter density perturbation and velocity divergence acquire the form 

4(t) = 6'^\t) + ^6^^\t) = 4^)(r) + I d'k,d^k2lCs{ki,k2;r)S^^!(.r)5i^^{r)SDiku - k), 

(3.1) 

(r) + I d3kid3k2/Ce(ki,k2;r)0i'^)(r)e£)(T)5,,(ki2-k), 

(3.2) 

= kj + kj + • • • and, following the procedure and 



where we introduced the notation k 



the results in 



we find that 



/C5(ki,k2;r) = /Cf(ki,k2) + 6/^L(ki,k2 



Mr) 



/Ce(ki,k2;T) 



C(ki,k2) + 6/^L(ki,k2 



' U2u2 
1^1 1^2 



(3.3) 
(3.4) 



where E{t) = Hl^rafi/ {^D+{a)). In the equations (^-(^) 



/Cf(ki,k2) 



5 2(ki •ks)^ ^ Iki •k2 {kl + kl) 
7^T 



+ 2- 



C(ki,k2) ^ ^ + - 



4 (ki • k2)2 



+ 2 



1 ki • k2 (A;? + 



^1^2 



are the standard one- loop Newtonian kernels for the expansion of the fields |38], and 

,5 - . 5 



/NL(kl,k2) 
/NL(kl,k2) 



(flNL - 1) + 0{1) 
(ONL - 1) + 0(1) 



(3.5) 
(3.6) 

(3.7) 
(3.8) 



are the terms arising from non-Gaussian initial statistics. Here onl parametrizes the pri- 



mordial NG level, according to the notations of Ref. [33|. The 0(1) terms in Eqs. ( p7i 
( |3.8| ) depend on the Fourier-space configuration. These terms represent general relativistic 
horizon-scale corrections |33| and have proven to be relevant in the description of the clus- 
tering of halos in the presence of non-Gaussian initial conditions as was recently discussed 
in [28|. Nevertheless, for our purposes here, we will deal with the primordial NG in the 
limit Ionl — 1| ^ 1, so, these terms become irrelevant and we will keep just the constant 
part of the non-Gaussian terms /^^ — 5(aNL — l)/3. ^ 

We assume that primordial NG is encoded in the curvature perturbations described through 
the Bardeen's gauge invariant potential <I> which, on sub-horizon scales reduces to the pe- 
culiar gravitational potential (/> up to a minus sign. The second order expansion of the 
matter density and velocity divergence can be expressed in terms of <I>, obtaining 

S^{t) = 7W(/£,a)$(k) 



^Notice also that we are neglecting the non-Gaussian contribution arising from the second-order evolution 
of perturbations during radiation dominance discussed in and since the matching at second-order 
into the radiation era will generate a term which then scale as a(r) cx r^, thus being subdominant on small 
scales w.r.t. the terms accounted here. 
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+ y"d3kid3k2/Cf(ki,k2)A^(A;i,a)4'^Al(A;2,a)cI)g,5D(ki2-k) (3.9) 

-^-^ = M{k,am\.) 

+ y"d3kid3k2/C^(ki,k2)A^(A:i,a)cI>£)A^(A;2,a)cI)g5B(ki2-k), (3.10) 

where the primordial gravitational potential ^ includes the linear part and the contribu- 
tions coming from primordial non-Gaussianities. The function M comes from the Poisson 
equation 

M{k,a) = ^^^D^{a), (3.11) 

with T{k) being the transfer function of matter fluctuations normalized such that T{k) 1 
when A; — )• 0. The expressions above account for a generic shape of the primordial NG, and 
also for a ACDM universe. Notice that, while for the primordial gravitational potential, 
the first lines of Eq. (|3.9| ) and ( |3.1C1| ) are exact, over the second lines applies the same 
approximations that led to eqs. (|3.1|) - (|3.2D . Using the expressions above, and recalling 
the definition ( |2.5[ ) in which we replace <^2 ^2/ /{t) in the ACDM case, the initial BS 
components can be written as 



(V?i(ki)(/?i(k2)99i(k3)) = e"^''(5D(ki +k2 + k3) x 

[2/Cf (ki ,k2)M\ki,a)MHk2, a)P$ (fci )P$ (^2) + cycl] 
+ e-^^M{ki,a)M{k2,a)M{k3,a) ($p(ki)cl>p(k2)<i>p(k3)) (3.12) 
(V?2(ki)(/?2(k2)¥?2(k3)) = e"^''(5z3(ki + ^2 + ks) x 

[2/C^ (ki ,k2)M\ki,a)M\k2, a)P$ (fci )P$ (^2) + cycl] 
+ e-^'iM{ki,a)M{k2,a)M{k^,a) {%{ki)%{k2)%{k^)) (3.13) 
(99i(ki)992(k2)v32(k3)) = 2e-3''5D(ki + k2 + k3) X 

[lCU'^iM)M\ki,a)M\k2,a)P^{ki)P^{k2) 
+ K.^{kiM)M\ki,a)M\k^,a)P^{ki)P^{k^) 
+ }Cf{k2,k3)M\k2,a)M\k3,a)P^{k2)P^{k3)] 
+ e-^^M{ki,a)M{k2,a)M{k3,a) {'^p{ki)%{k2)%{k3)) {3.U) 
(99i(ki)99i(k2)v32(k3)) = 2e-3''dD(ki +k2 + k3) X 

[lC^{ki,k2)M\kua)M^{k2,a)P^{ki)P^{k2) 
+ fC^{ki,k3)M\ki,a)M\k3,a)P^{ki)P^{k3) 
+ }C^{k2,k3)M\k2,a)M\k3,a)P^{k2)P^{k3)] 
+ e-^^ M{ki,a)M{k2,a)Mik3, a) {'^p{ki)'^p{k2)%ik3)) (3.15) 

where P^{k) is the gravitational potential power-spectrum, see Eq. ( ^I20|) below. The 
expressions above constitute the initial conditions for the bispectrum in the system ( 2.11] ) 
and must be evaluated at the initial time, for which rj = loga/aj„, = 0. Notice that the 
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expressions above are separated in a perturbative Newtonian term, which is the first term 
on the r.h.s of each equation, and in a primordial part, which is the term including the 
three-point correlator of the primordial gravitational potential. 



3.1 Primordial NG: Shapes and Running 

The functional form of the bispectrum of the primordial gravitational potential entering in 

Eqs. (ra)-(ra)- 

($(ki)$(k2)$(k3)) = 5D(ki + k2 + k3)B^{h,k2, ks) (3.16) 

depends on the details of the mechanism that generated the primordial fluctuations. In 
the following we will analyze three representative phenomenological models associated to 
different mechanisms generating the primordial NG. Note that there exist other possible 
configurations (see for instance |^l[) which correspond to more general deviations from the 
standard slow-roll inflation and which are not well described by the models treated here. 

3.1.1 Local Shape 

In the local model, the NG for Bardeen's gauge invariant primordial gravitational potential 
$ is generated by a quadratic expansion which is local in real space 

a>(x) = $l(x) + /i^r'($^(x) - ($2)) , (3.17) 

where is the linear Gaussian part of the potential. The dimensionless constant /nl 
defines the NG strength in this model. One must also notice that, since $ and evolve 
proportionally to g{a) = D^{a)/a, the Eq. ( 3.17] ), and consequently, the definition of 



/nL) depends on the choice of the redshift at which this equation is extrapolated. The 
relation between two extrapolation choices a and b is expressed through /^^ = = 
zjj)/g[z = Za)f^j^. Two common conventions in the literature are the LSS convention in 
which Eq. (3.17) is extrapolated to z = 0, and the CMB convention which extrapolates to 



z = oo. They are related through f^f^ = g{z = oo)/ g{0)f^f^^ where g{z = oo)/g{0) ~ 
1.3064. In this paper, we assume the CMB convention, which implies that $ is evaluated 
deep in the matter era. From equation ( |3.17 ) one can easily derive the expression for the 
local bispectrum 

B^{ki,k2, ks) = fhT'F'°'"\ki,k2,k3), (3.18) 

where 

A:2, ks) = 2 [P^{ki)P^{k2) + P^{ki)P^{ks) + P^{k2)P^{k3)] (3.19) 

, 1 1 1 

2A| • ( . + . + 



,(/fciA;2)3-("»-i) (fciA:3)3-("=-i) (A:2fc3)3~("=-i) , 
and P^{k) is the gravitational potential power spectrum 

($(ki)«>(k2)) = 5^{ki2)P^{ki) = 53(ki2)A$ • A:~3+(n.-i) (3 
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with scalar amplitude A$ and tilt Ug- In slow-roll single scalar field inflation, /jIJl^' is 
suppressed by the slow-roll parameters and consequently the primordial NG is unmeasur- 



ably small [13, 14 1. However, in the presence of light fields other than the inflaton, as for 
instance in the curvaton model or in multi-field models, inflation can produce large NG of 
the local type with I/^l^'I ~ 0(10^) |15|. The strongest current limits on the strength of 



local NG set the f^{] parameter to be in the range —4 < < 80 at 95% confidence level 



3.1.2 Equilateral Shape 

We can also consider inflationary models with higher order derivative operators and non 
canonical kinetic terms, such as the DBl model |44], for instance. In this case, as was 
discussed in the primordial bispectrum is well described by the equilateral template 
Bis,{ki,k2,k3) = /^'[F^^(A;i,fc2, fea), where 

F''i{ki,k2,k3) = 6 (-P^{k^)P^{k2) + (2 perm.) - 2 {P$(A:i)P$(fe2)P*(A;3)}'/' 



+ ( P^iki) {P$(fe2)}' {P'S'ih)}' ) + (5 perm.) 



6^1 • f -vnr^Lf;^^ - (2 perm. 



(fcl^2)3-("-l) ^ " ^k,k2k3f-¥-^'^) 
1 



2-|(n,-l) 3-K-l) 



+ (5 perm.) . (3.21) 



"1 "-2 "'3 

The primordial bispectrum is maximized by configurations with modes of similar momen- 
tum scales ki k2 ^ k^, i.e. for equilateral configurations. At present /^^ is constrained 

feq 
'NL 



to be -125 < /^l < 435 (95% C.L.) |3 



3.1.3 Folded Shape 

A third template is the so called folded or flattened model. This model is related to the non- 
Gaussianities generated by deviations coming from the choice of a non-standard adiabatic 



Bunch-Davies vacuum state as initial state |46, 47]. The associated primordial bispec- 



trum assumes a complicated functional form in terms of the momenta /cj. Nevertheless, 
in Ref. it is proposed a factorized template which reflects very well the main features 
of the primordial NG associated to the choice of a non standard vacuum. The proposed 
factorized template is B^{ki, k2, ks) = f^^F^°^'{ki, k2, k^), with 

F'°Hkuk2,k3) = 6 [p^{k^)P^{k2) + (2 perm.) + 3{P<i,(A;i)P$(fc2)P$(A;3)}2/3 

1 /3 

P$(fci) {P$(A:2)}' {P^{k3)f) + (5 perm.; 



6A|- I 77-17^^7-71 + (2 perm.) + 



(A:iA:2)3-K-i) " {kik2k3f-¥^'-'^ 
1 



l-lK-l) 2-§(n,-l) 3-(n.-l) 
"-l '^2 '''3 



(5 perm.) . (3.22) 
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kl 

(1) Squeezed (2) Equilateral (3) Folded 

Figure 4: Momentum triangle configurations which maximize the templates functions F{ki, fc2, ^a)- 
(1) The squeezed configuration fcs ^ fci « ^2 maximize the template function F^°'^^^{ki, k2, k^) 
produced in local models of inflation. (2) The equilateral configuration fci w /c2 ~ ^3 maximize 
F°'^(ki^k2,k'i) which arises on models with higher-order derivatives and non canonical terms. (3) 
The "folded" or "flattened" configuration k2 ~ k^ ^ fci/2 maximize the F^°' (fci, fc2, ^a) which is 
related to the choice of a non-standard vacuum other than the Bunch-Davies vacuum. 

The primordial bispectrum for this model is maximized by configurations with modes 
obeying the momentum configuration A:2 ~ ^3 ~ ^i/2. The triangle configurations relevant 
for each model of NG are represented in figure Q 

3.1.4 Running /nl 

So far, we have considered a scale independent /nl parameter regardless the shape of 
the primordial bispectrum. However, recently, several inflationary scenarios have been 
considered which naturally allow for NG with a generic scale dependence on /nl- More 



specifically, it was shown in Refs. EH] and |49| that single scalar field models of inflation 



with variable speed of sound Cg, predict in general a primordial bispectrum of the equilateral 
type with a scale dependent /nl- In the following, for non-Gaussianities of the equilateral 
type, we adopt the functional form for the scale dependence proposed in []50| : 

hh{ki,k2,k^) = /nl,p {~^~^^~~^ ' 

where fT<sh,p is the non linear parameter evaluated at some pivot scale kp and k is a free 
parameter (related to the speed of sound) that quantifies the running on the scale and 
which is constant at least between CMB and cluster scales [^]. To be consistent with 
observational constraints, the scale dependence must be very mild, then \k\ <C 1. Also, 
in an interesting discussion regarding the scale dependence of the biasing parameters for 
the equilateral model, in Ref. |^5|, the authors found the constraint \k\ < 0.3 in order to 
guarantee the convergence of the integrals quantifying the scale dependence of the biasing. 
In this paper we thus consider values of the running parameter within \k\ < 0.3. Moreover, 
we use negatives values of k since in this case the NG is enhanced for scales beyond the 
pivot scale Another proposal for the functional dependence on the running /nl was 



discussed recently in [^| 



/NL(A;i,fe2,A:3) = /NL,pf^^4^1 ' (3-24) 
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which replaces the arithmetic mean of the momenta in equation ( |3.23| ) with their geomet- 
ric mean. Notice that the two functions coincide for equilateral configurations. While the 
last proposal has the virtue of being separable, which simplifies significantly the numerical 



implementation of the CMB estimator, the equation (3.23) provides a more accurate de- 
scription of the primordial bispectrum in DBI models. The functional dependence (3.24) 
seems to be better suited for NG of the local type. It was shown recently that non- 
Gaussianities of the local type exhibit significant scale dependence of the amplitude of the 
primordial bispectrum |5^]. The form of this amplitude depends on the specific details 
of the underlying inflationary model. For example in the case of standard single-field mod- 
els of inflation it coincides with the functional form ( |3.24 ) in the equilateral limit. Away 
from this limit, it is not possible to factor out the geometric mean of the momenta in 
/nl but, in principle, the essentials of the local model can be roughly captured with the 
expression ( 3.24| ). 



4. Results 



For our analysis, we assume cosmological parameters corresponding to a A-CDM model in 



agreement with the WMAP 5-year data [42|. The cosmological parameters for this model 
are summarized in Table 1. As discussed before, we use second order perturbation theory 
from z = oo to z = Zin, and the TRG from Zin to z = 0. To get the numerical results 
presented in the following, we set Zm = 50, and used the CAMB code [43| to derive the 
linear power-spectrum. In order to reproduce the cosmological model in table 1, we must 
flx the scalar amplitude to be = 2.48 x 10~^ in the notations of Ref. [42]. With these 
parameters we are able to compare our results for the matter power spectrum with the 
results of N-body simulation of cosmological structure formation with non-Gaussian initial 

M for N-body 



conditions of the local type from Refs. [|5J] and |Q (see also ||5^, |2^, [ 
simulations with local-type non-Gaussian initial conditions). 

We also compute the reduced bispectrum for matter density perturbations. Here we 
point out some general features about the effect of the primordial non-Gaussianities on the 
bispectrum based on our results. 



Table 1: A-CDM cosmological parameters assumed on this paper. 



Model 


h 


0-8 


Us 








A-CDM WMAP5 


0.701 


0.817 


0.96 


0.279 


0.721 


0.0462 



4.1 The Power Spectrum 

Here, we present the results of the power spectrum of the matter density perturbations 
resulting from the integration of the system ( 2.11| ) on the scales of interest for next gener- 
ation galaxy surveys. In order to evaluate the impact of primordial non-Gaussianities on 
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the power-spectrum we plot the ratio of the power spectra for non-Gaussian and Gaussian 
initial conditions for different values of the /nl parameter. It must be noticed that the 
power spectra for both Gaussian and non-Gaussian initial conditions were calculated using 
the TRG method. 



4.1.1 Local Model 



For the local non-Gaussian model, the result of this evaluation is shown in figure |. For 
this type of NG we can compare our findings with the results of N-body simulations in p^ . 
In the figures we also show the result from one-loop perturbation theory, see for instance 
|35]. We notice that up to ~ 0.2/i/Mpc, the one-loop and the TRG describes the data 
very accurately while, for larger values of k we begin to have departures from the data. 
The difference between one-loop and TRG prediction becomes noticeable for large /nl and 
for higher wavenumbers and lower redshifts as evidenced in the right panel of Figure |5|, 
which is the result of the evaluation at z = 0. From the figures we can also infer that 
the TRG approach allows a less suppressed non-linear growth of the power spectrum in 
the region k > 0.25/i/Mpc. This tendency seems to be general for the different models 
of primordial NG. The oscillatory behaviour exhibited in the TRG plots are the result of 
numerical errors in the numerical integration of the system p. 11 . 





TRG 








1-loop 








NBody 


Local 








z=1 








fnfcSOO 










^ fnl=250..- 


' ■ ' fnl=80 

1 ■ ■ 




k (h/Mpc) 



k (h/Mpc) 



Figure 5: Ratio of the non-Gaussian to Gaussian power spectrum for several values of /nl in 
the local model. The dots correspond to the data from the N-body simulations of |5^. The red 
(continuous) line is the TRG result of this paper and the blue (dashed) line is the one-loop result. 



4.1.2 Equilateral and Folded Shapes 

Here, we show the results of the power spectra ratio for equilateral and folded non-Gaussian 
shapes. The result of this evaluation is shown in Figures ^. Since at present there are no 
simulations using this type of primordial NG, the results displayed here have a predictive 
character. Again, as we observed in the local case, when we compare with the one-loop 
calculation the TRG predictions allow for less suppressed effects of non-Gaussianities, The 
difference is more evident for scales smaller than k ^ 0.2/i/Mpc and for lower redshifts. 
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Figure 6: Ratio of the non-Gaussian to Gaussian power spectrum for several values of /nl in the 
equilateral (top panels) and folded (bottom panels) models. The red (continuous) lines are the 
TRG result of this paper and the blue (dashed) lines are the one-loop result. 



4.1.3 Running /nl? equilateral model 

In Figure |^ we plot the power spectrum ratio for the equilateral model with scale de- 
pendent according to Eq. ( 3.25 ). We choose the pivot scale kp = 0.04/i/Mpc, and 
K = 0, 0.25, —0.25. We see that for positive k we have a suppression in the power spectrum, 
while, for negative k we can get a significant enhancement of the power spectrum. No- 
tice that the scale dependence of /^^ can make the effects of primordial NG in the power 
spectrum of the same order of the local case. 



5. The Reduced Bispectrum 

As we will be mainly interested in the study of the dependence of the bispectrum on the 
shape of the momentum triangle, it is convenient to focus the analysis on the reduced 
bispectrum which is a useful quantity defined as^ 

Q[ki,k2, k-i) = (5.1) 
S(A;i, fc2, fcs) 



^For brevity, in the following we omit the time dependence and the subscripts in the correlators since 
we will deal only with the matter density component. 
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■ TRG 

■ 1-loop 



Running-Equilateral 
fnl=250, z=0 




k (h/Mpc) 



Figure 7: Ratio of tlie non-Gaussian to Gaussian power spectrum for /nl = 250 in tlie equilateral 
model for different values of k. The red (continuous) lines are the TRG result of this paper and the 
blue (dashed) lines are the one-loop result. 



where T,{ki, k2, h) = P{ki)P{k2)+P{ki)P{k:i)+P{k2)P{k:i). The greatest challenge in the 
interpretation of galaxy clustering data from any surveys is galaxy bias. As thoroughlly 
discussed in Ref. [26|, the reduced bispectrum provides an excellent determination of 
linear and non-linear bias parameters of intermediate and high-redhift galaxies, when all 
measurable triangle configurations down to mildly non-linear scales are included. The 
reduced bispectrum is also a powerful probe of primordial NG. The planned galaxy surveys 
at redshift z > 2 should yield constraints on NG that are comparable to, or even better 
than, those from CMB experiments. The precise computation of the reduced bispectrum 
is therefore of primary importance. 

Due to statistical homogeneity and isotropy, the bispectrum and the reduced bispec- 
trum depends on time, on the magnitude of two of the momenta, for instance, ki and /c2, 
and the angle between them cos 9 = ki ■ k2. Figure ^ shows the geometric parameters in 
the three configurations that we consider in this paper. For gaussian initial conditions and 
an equilateral configuration, Q{k, k, k) ~ 0.57 at the tree-level in perturbation theory and 
is independent of scales. On the other hand, Q{k,k,k) exibits a clear scale dependence 
when some NG is present. The same is true when one departs from the equilateral configu- 
ration which reflects the anisotropy on the growth of structures dictated by the non linear 
Eulerian dynamics. As the NG appears in our RG equations as post-Newtonian term and 
therefore suppressed at small scales and late times compared to the leading Newtonian 
terms, it is expected that the NG will be more easily testable at high redshifts and on 
large scales. This is particularly true for the local shape as the reduced bispectrum (at the 
tree-level) which is directly proportional to scales like 1/Ai{k,a). 

In our analysis we fix the ratio k2/ki = 2 and evaluate for different values of ki, /nl 
and for different redshifts. As an extreme case, we also consider the ratio k2/ki = 1 to have 
some indications of the maximization of the non-Gaussian effects in the squeezed limit. 
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5.1 Local model 

The results of our evaluations for the local model of primordial NG can be seen in Figure 
Hj. From the upper panels in Figure ^, which corresponds to A; = 0.01 /i/Mpc, we see that 
Q is much more sensitive to the effect of the primordial NG for small ki. In this case 
we see that the effect of non-Gaussianities on Q is almost independent of 6 for z = \ 
and z = 0. Nevertheless, for higher redshifts we can see clearly a maximization in the 
effect of the non Gaussianities in the squeezed limit which corresponds to ^ — )• tt. The 
effect is dramatically maximized when we consider the squeezed limit of the configuration 
^2/^1 = 1 as shown by the thin lines of this Figure. In this configuration, the condition 
— 7> TT is equivalent to k'i/ki <C 1. We can understand the behaviour at this limit by 
studying the first terms in the perturbative expansion of Q. The first term containing 
the non-Gaussian contributions is A:2, /cs) = B^^'>^^'^{ki,k2,k?,)/TP{ki,k2,k:i). 

Going to the squeezed limit k^/ki <^ 1, this non-Gaussian contribution scales roughly as 
q(o),ng^^^^ /ci, A;3) ~ /nl/^3- This is why in the limit — > vr, the non-Gaussian component 
is huge, actually, it is divergent. For this reason we should not consider the evaluation at 
this extreme configuration as a precise prediction but as an indication of the behaviour 
at this limit. The exact scaling of Q^^^'^^ depends also on the precise details of the 
transfer function which could alleviate up to some extent the divergent behaviour in this 
configuration for higher values of ki as suggested by the results in the lower panels of 
Figure ^ On the other hand, the tree level Gaussian component of the reduced bispectrum 
Q(0)(A;i,A:2,A:3) = S*^"^(A;i, ^2, A;3)/S0(A;i, fea, ^3), where 5*''^"(A;i, fes, ^3) is the Newtonian 
piece of the bispectrum, scales like Q^^\ki,ki,k3) ~ k^/ki in the squeezed limit and 
consequently it goes to zero for 9 ^ tt. This tendency can be seen in the black-thin lines 
in the upper panels of Figure ^. Finally, the impact of non Gaussianities for higher values 
of ki is much smaller and is almost independent of /nl as can be seen in the bottom panel 
of Figure ^ which corresponds to A: = 0.1 h/Mpc. 

Figure ^ shows the comparison of the TRG evaluation with the perturbation theory 
(PT) tree level reduced bispectrum Q^^^ for several values the /nl at redshift z = 0. These 
plots are meant to help understanding the relevance of non linear effects at different scales 
and redshifts. It is noticeable that there is a remarkable change in the impact of non linear 
effects depending on the scales we are considering. For instance, in the left panel of figure 
m, we plot for ki = 0.01/i/Mpc. At this scale, the /nl dependence is strong, as we can 
see from the large separation in the solid lines which depends of the strength of primordial 
non-Gaussian effects. On the other hand, the effects from non-linearities quantified by the 
separation between dashed lines and the solid lines is rather mild. In the right panel, we 
plot the same situation for ki = 0.1/i/Mpc. At this scale, we observe that the situation 
is quite the opposite, the non-linear effects are stronger, while, the dependence on the 
primordial non-Gaussianities is much weaker. 



5.2 Equilateral model 

The equilateral model exhibits a different pattern compared to the local model. This 
behaviour can be seen more clearly in the upper panel of figure 10. There, we see that Q is 
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Figure 8: The reduced bispectrum for the local model of NG with fixed k\ and ratio ^2/^1 = 2 
(thick lines) and fc2/fci — 1 (thin lines). The left panels shows the evaluation at z = 1 while the right 
panels does it for z = 0. In each case we plot Q for non-Gaussian initial conditions with /nl = 250 
(green-dotted) and /nl — 500 (red-dashed). From top to bottom we plot Q for fc = 0.01, 0.05, and 
0.1/i/Mpc respectively. The thin red-dashed and green-dotted lines in the upper panel show the 
divergent behaviour due to the non- Gaussian component near = tt in the squeezed limit. The 
black-continuous thin line on the same panels represent the Gaussian component of the reduced 
bispectrum which approaches to zero as tt in this configuration. 

enhanced by non-Gaussianities in a way that it is highly dependent on Q. The maximum of 
the effect occurs at ~ 27r/3 which corresponds to equilateral configurations in this model 
and are responsible for the bump in this Figure near to O.Svr. Again, as in the local 



-19- 




Figure 9: Comparison of tlie tree level Q^^^ and the TRG result for the reduced bispectrum. 
Continuous-thick lines are the TRG result and dashed-thick lines represent the tree level evaluation 
at the configuration k2/ki = 2. Both evaluations are done at z = for different values of /nl and 
for fci = O.Olh/Mpc (left), and fci — O.lh/Mpc (right). Continuous and dashed thin lines represent 
the evaluation for k2/ki — 1. 



model, the effect of non-Gaussianities is greatly enhanced in the configuration k2/ki = 1 
while for higher ki the reduced bispectrum is much less sensitive to the effect of primordial 
NG regardless the value of the /nl parameter. 



The result of the evaluation of the running /nl equilateral case is shown in figure 11 
In the configuration kijkx = 2, we evaluate for k = 0,0.25,-0.25. For ki = 0.01/i/Mpc 
and z = 1 we find a remarkable enhancement (suppression) for positive (negative) values 
of K. For ki = 0.1 /i/Mpc and z = the effect is rather unnoticeable. 

5.3 Folded model 

Contrary to the equilateral model, in the folded model, Q is enhanced by the non-Gaussian 
contributions coming from the correlations among modes with almost coUinear wave- vectors 
{9 = 0, vr). This effect is maximized in the limit 6 ^ ir (as in the local model) which 
corresponds to a folded configuration of the momentum-space triangle when k2 ~ fes ~ A;i/2 



(Figure Q). This effect can be seen in figure 12. As in the previus cases, the effect is 



maximized for small ki and for high redshift. For higher ki and lower redshift as in the 



bottom panel of figure 12 we see that the effect of this model of non-Gaussianities is tiny. 



regardless the value of /nl- 



6. Conclusions and Discussion 

In this work we have studied several aspects related to the effects of primordial non- 
Gaussianities in gravitational clustering using the recently developed TRG method. We 
have discussed different shapes of primordial NG and evaluated their impact in the non- 
linear evolution of the power-spectrum and (reduced) bispectrum of matter density per- 
turbations. 
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Figure 10: The same as figure 8 for the equilateral model of primordial NG. 



For the power-spectrum with some primordial NG present, we have compared our re- 
sults to the ones of standard one-loop perturbation theory and N-body data. In the range 
0.1 < kMpc/h < 0.2, the one-loop and the TRG method reproduce pretty well the data 
from the N-body simulations of Ref. ||5^ while, for k > 0.2 /i/Mpc, the one-loop prediction 
begins to show some inaccuracies which are within the ^ 2 % level. The non-linear growth 
of the power-spectrum is more suppressed in the one-loop calculation than in the TRG 
calculation. The differences between the two approaches becomes more evident on smaller 
scales and lower redshifts, where the TRG offers a noticeable improvement with respect 
to the one-loop but still with a deviation from the N-body data. We expect that such 
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Figure 11: Reduced bispectrum for running /nl with /nl,p = 250 and k — (solid, black line), 
K = 0.25 (dashed, red line) and k — —0.25 (dot-dashed, blue line). 



deviations can be reduced if we take into account the "vertex corrections" including the 
trispectrum in the system of equations as we discussed in the Introduction. On general 
grounds, we expect that the effects arising from vertex corrections will become relevant for 
higher wavenumbers and that they will help to extend the range of validity of wavenumbers 
to larger values. We have also computed the reduced bispectrum for various shapes of NG. 



As stressed in Ref. [26|, the reduced bispectrum is a powerful probe of NG. There the au- 
thors performed a Fisher matrix analysis to study the smallest value of the paramater /nl 
measurable in high-redshift galaxy survays for local and equilateral shapes, after marginal- 
izing over the bias parameters. At the tree-level the signal-to-noise ratio (S/N)'^ scales 
like (the sum over momenta forming a triangle of) B'^{ki, k2, k^)/ P{ki)P{k2)P{k3). For a 
local shape of NG, A;i < A;2 ~ ^3, B'^{ki,k2,k^) / P{ki)P{k2)P{k^) ~ P{ki). Non-linearities 
enter in the determination of the maximum wavenumber A;niax at which one computes 



(S/N) ~ /cmax- For the equilateral shape, fci ~ ^2 ~ /^s, and again at the tree level, the 



signal-to-noise ratio scales like (S/N) ~ /^l ^max -P"'^''^(^max) /NL^max- Going beyond 
the linear order through the TRG method one should compute the full Fisher matrix and 
therefore (5-B(fei, A;2, A;3)/9/nl)^- In this case the non-linear effects cannot be neglected. 
We leave this and other issues, such as the inclusion of the trispectrum, for future work. 
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Figure 12: The same as figure 8 for the folded model of primordial NG. 
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